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I. INTRODUCTION 

When quantum particles in a periodic potential of pe- 
riod d are subject to a weak additional constant force F 
they do not uniformly accelerate like free particles, but 
instead undergo Bloch oscillations p] at an angular fre- 
quency given by: 

uj b = Fd/h . (1) 

Bloch oscillations (BOs) of cold atoms in optical lattices 
were first observed in 1996 by uniformly accelerating the 
lattice 0: in a frame co-moving with the lattice the 
atoms experience a constant force. At about the same 
time, the accelerating lattice method was used to observe 
Wannier-Stark ladders [3] , which are a different aspect of 
the same "tilted lattice" physics. The method has subse- 
quently been employed to realize beam splitters for atom 
optics capable of large momentum transfers, see, e.g. [4]. 

In gravity-driven BOs the lattice is held fixed in space 
but oriented vertically so that gravity provides the force 
F g = mg on atoms of mass m. From Eq. (JT|, a mea- 
surement of wb corresponds to a measurement of the 
applied force F if we know d/h. This Bloch oscilla- 
tor may be viewed as an interferometer in momentum 
space [5] and has been experimentally demonstrated by 
a number of groups [BHS]. For example, the experiment 
[9] used gravity-driven BOs of strontium atoms to mea- 
sure the local acceleration due to gravity at the level of 
Ag/g = 5 x 10~ 6 . Like any interferometer, long coher- 
ence times are crucial for precision measurements and 
in [5] the BOs were coherent over 7 s, corresponding to 
« 4000 oscillations. This remarkable degree of coherence 
was greatly facilitated by the choice of strontium atoms, 
which have very weak s-wave scattering, and thus dy- 
namical instabilities normally associated with superflow 
in lattices jTUl [TT] were highly suppressed. Variations on 
this scheme that improve the visibility of the BOs, includ- 
ing frequency [12] and amplitude [Ml [14] modulation of 
the lattice, have allowed for the measurement of gravity 



at the level of Ag/g = 10 9 . In these latest experiments 
the BOs were coherent for over 20 seconds. 

The experiments referred to above all involve destruc- 
tive measurements of the BOs due to the nature of the 
imaging process of the atoms, whether it be in situ or 
by a time-of-flight technique after the lattice has been 
switched off [13]. Therefore, a precision measurement of 
wb by the above methods requires that the experiment 
be re-run many times, each run being for a slightly differ- 
ent hold time, so that the oscillations can be accurately 
mapped out. This not only takes a long time, but also 
requires that the initial conditions be recreated as faith- 
fully as possible for each run. 

In [T5] we proposed a scheme for continuous (i.e. non- 
destructive) measurements of BOs based upon placing 
the atoms inside a Fabry-Perot optical resonator which 
would allow for an estimate of uj-q from the data acquired 
over a single run. A related scheme has also been inde- 
pendently proposed for ring cavities |16j . The periodic 
potential is now provided by the standing wave of light 
which forms inside the cavity when it is pumped by a 
laser. Orienting the cavity vertically, the atoms execute 
BOs along the cavity axis as depicted in Fig. [T] The en- 
hanced atom-light coupling inside a high-Q cavity means 
that the oscillating atoms imprint a detectable periodic 
modulation on both the phase and amplitude of the light 
which can be seen either in transmission or reflection. 
Thus, the measurement is performed upon the light leak- 
ing out of the cavity rather than directly upon the atoms. 

The strong atom-light coupling that can be realized 
in cavity-QED stands in contrast to the case of optical 
lattices in free space where the atoms exert only a tiny 
backaction upon the light. The optical dipole interaction 
between a single cavity photon and a single atom is char- 
acterized by the Rabi frequency g — (fJ,/H) \fhuj c / \sqV), 
where w c and V are the frequency and volume of the rel- 
evant cavity mode and fj, is the atomic transition dipole 
moment. Defining the cooperativity C = <7q/(27k), 
where 2j is the spontaneous emission rate of the atom 
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in free space and 2k is the energy damping rate of the 
cavity, 1/C is the number of atoms required to strongly 
perturb the light field. The normal mode splitting that 
results from strong coupling has been directly observed 
in a number of cold atom optical cavity experiments 
[T7H20] . In the experiment [20] , which was performed 
with a Bose- Einstein condensate, the cooperativity was 
C = 145. Even more pertinently, these systems have 
been used to detect the presence of single atoms [2"TH2"3"] , 
as well as to follow their dynamics in real time (24) [25] . 
The collective dynamics of ultracold atomic gases have 
also been tracked using cavities . The key experi- 

mental steps necessary for the continuous monitoring of 
BOs in a cavity have, therefore, already been demon- 
strated. 

The rub with any continuous measurement scheme is 
measurement backaction. In cavities this backaction typi- 
cally takes the form of cavity photon number fluctuations 
which lead to random force fluctuations on the atoms, as 
is evident in the erratic nature of the single atom trajec- 
tories seen in the experiments [24, 25; referred to above. 
In the many atom context, quantum measurement back- 
action generally manifests itself in a heating of the atom 
cloud (although under some circumstances it can lead 
to cooling 129]). In the cavity-optomechanical regime 
(where the collective motion can be modelled as a har- 
monic oscillator of angular frequency ui) the heating rate 
is expected to be R = (x zp {/h) 2 Sjrjr(—uj) [3U], where 
x zp { is the zero-point fluctuation and Sj^jr is the spec- 
tral density of the force fluctuations (which is directly 
proportional to the cavity photon number fluctuations). 
This heating rate is in agreement with observations when 
convolved with technical fluctuations [3"T] . 

In the system considered in this paper (see Fig. [T]) , we 
can divide the backaction into two types. One type comes 
from the fact that the atoms sit in an optical lattice whose 
depth is periodically modulated in time at the frequency 
wb due to the effect of BOs. This backaction is a classical 
effect in the sense that it occurs even when the light field 
is treated classically (no photons) . The nonlinearity that 
arises from this backaction can lead to swallowtail loops 
in the atomic band structure |321 134) that mimics the 
effects of direct atom-atom interactions 35-38]. These 
loops are the counterpart in the atomic wave function 
of optical bistability in the light |2l)H2gl 153] . The sec- 
ond type of backaction arises only when the fluctuations 
due to the discrete photon nature of the light field are 
taken into account and is related to the heating effect 
mentioned above. The characteristic frequency of these 
latter fluctuations is k which is much larger than wb- 

The first type of backaction was analyzed in our pre- 
vious paper [15] where our main aim was to show that, 
despite the self-generated time modulation of the intra- 
cavity optical lattice, the Bloch acceleration theorem still 
applies and the BO frequency is not modified (although 
harmonics can be generated) . This latter result is clearly 
very important if the cavity BO method is to be used for 
precision measurements and may be viewed as a conse- 




FIG. 1. (Color online) Schematic of the precision measure- 
ment proposal in [T5]. A dilute cloud of cold atoms under- 
goes BOs in the combined intracavity lattice potential and 
the acceleration due to gravity. The transmitted light field's 
intensity and phase are modulated at the Bloch frequency. 
A in-situ precise measurement of the Bloch frequency (and 
hence the force) can be performed by detecting the transmit- 
ted light. 

quence of the fact that the formula ([I]) does not depend 
on the depth of the lattice, only its spatial period. An es- 
timate of the effects of the second type of backaction was 
also given in our previous paper, but this estimate was 
obtained under the assumption that the photon number 
fluctuations were purely due to the photon shot noise 
found in a coherent state of light. This ignores the cor- 
relations that build up between the atoms and the light 
inside the cavity and our main aim in this paper is to 
solve the dynamics of the coupled photon and atom fluc- 
tuations systematically from first principles and thereby 
capture these correlations. This will allow us to prop- 
erly determine the sensitivity of the measurement of the 
Bloch frequency to quantum fluctuations. 

The plan of this paper is as follows: in Section[n]we in- 
troduce the physical system, the associated hamiltonian, 
and the equations of motion. We then review in Sections 
|III| and |IV| the meanfield approximation and the associ- 
ated numerical results which were the focus of our previ- 
ous paper |15j . before introducing in Section [V] the main 
model to be treated in this paper which adds quantum 
fluctuations. This is an elaboration of the linearization 
approach presented in, e.g. [39H4T] . to include a time- 
dependent meanfield component (due to the BOs). The 
fluctuations correspond to quasiparticles (excitations out 
of the meanfield), and their spectrum is analyzed in Sec- 
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tion VI and then used to help interpret the numerical 



results for the quantum dynamics presented in Section 
|VII| We also develop a simple rate equation picture, valid 
in the weak coupling regime, to help us understand the 
rate of quasiparticle excitation. Following this we change 
gears slightly and apply the above results to investigate 
how quantum fluctuations affect a precision measurement 
of wb by calculating the signal-to- noise ratio (SNR). We 
present the theory of lying behind these calculations in 
Section |VIII| and in Section |IX| we examine the results, 
paying particular attention to whether or not there is 
an optimal value for the atom-light coupling parameter 
/3 = NUo/k. We also present results illustrating the de- 
pendence of the SNR on other system parameters such 
as the number of atoms and the intracavity lattice depth. 
We summarize our results and provide some perspective 
in Section [X] We have also provided two appendices that 
give details omitted from the main text, including those 
of our numerical approach. 



II. HAMILTONIAN AND EQUATIONS OF 
MOTION 

Our system consists of a gas of N bosonic atoms inside 
a vertically oriented Fabry-Perot optical cavity. A single 
cavity mode of frequency uj c is externally pumped by co- 
herent laser light with a frequency w p that is detuned 
from the atomic resonance and the cavity resonance fre- 
quency. This sets up a standing wave mode along the 
cavity axis of the form cos(fc c z), where k c = u> c /c. The 
relevant frequency relations are characterized by the two 
detunings 



A c = uj p - u> c 

A a = 0J n - LU n 



(2a) 
(2b) 



where uj a is the atomic transition frequency. In the dis- 
persive regime, the occupation of the excited atomic state 
is vanishingly small and it can be adiabatically elim- 
inated. A onc-dimcnsional hamiltonian for the atom- 
cavity system in the dispersive regime can then be writ- 
ten as HI EH 



H 



hA c a* a + ihrj (a) - 

h 2 d 2 



dz * 



2M dz 2 



+ HUoa'a cos 2 (k c z) — Fz 



(3) 



where &(z, £) and a(£) are the field operators for the 
atoms and the cavity photons which obey the equal time 
bosonic commutation relations [9(x, £), W{pd , £)] = 5(x— 
x'), and [&(£),&*(£)] = 1, respectively. The single atom 
dispersive light shift has been denoted by Uq = g 2 /A a . 

The hamiltonian has been written in a frame rotat- 
ing with the pump laser frequency cj p , and this leads to 
the appearance of the two detunings. The first term is 
just the free evolution of the cavity mode. The second 



term represents the laser coherently pumping the cavity 
at rate ry, and the third term describes the atomic part of 
the hamiltonian. The first two terms of the atomic part 
represent the kinetic energy and a light induced poten- 
tial energy. This latter term can either be understood 
as the atom moving in a periodic potential with average 
amplitude HUo^a) or, if combined with the first term 
in the hamiltonian, as a shift in the resonance frequency 
of the cavity due to the coupling between the atom and 
the field. The third term in the atomic part provides 
the external force that drives the BOs. We assume this 
force arises from the vertical orientation (z increases in 
the downward direction) of the cavity and is given by 
F = Mg. 

Note that we have not included direct atom-atom in- 
teractions. From a practical point of view this is because 
these are easily dwarfed by the atom-cavity interactions, 
or can be tuned away using a Feshbach resonance |46) . 
However, from a more pedagogical point of view, we want 
to focus here upon the nonlinearity which arises from 
atom-cavity interactions and avoid any confusion with 
the effects of direct atom-atom interactions which are in 
any case already well studied in quantum gases. The 
interactions between the atoms and the common light 
mode provided by the cavity gives rise to a nonlinearity 
that is in some ways analogous to the direct interactions 
but in other ways differs and can lead to novel behavior 



Natural units for the length and energy in cavity-QED 
are given by l/k c and the recoil energy En = h 2 k 2 /(2M ), 
respectively. From here on we scale all lengths by l/k c 
and consequently define x = k c z. We scale frequencies by 
the recoil frequency wr. = E^/h and time by 1/wr and 
retain the same symbols for the scaled variables. The 
Heisenberg-Langevin equations of motion for the light 
and atomic field operators in the scaled variables are S3] 



da 

! — 

dt 



dx ^(x,t)^(x 1 t)U Q cos 2 (x) - in a 
+ irj + iV2ni(t) (4a) 



dx 2 



+ JToa^acos 2 ^) — fx 



(4b) 



where / = F/{hk c ujYi) — wb/(7twr,) is the dimension- 
less form for the external force. The operator £(£) is the 
Langevin term and is assumed to be Gaussian white noise 
with the only non-zero correlation being 



(£(#(£')) =*(*- O- 



(4c) 



Mathematically, the Langevin noise terms are neces- 
sary in order to preserve the commutation relation 
[a(£),a^(£)] = 1 in an open system. Physically, their ori- 
gin is vacuum fluctuations of the electromagnetic field 
that are transmitted into the cavity via the mirrors and 
they thus only appear in the equations for the light field. 
Nevertheless, the noise is conveyed to the atomic dynam- 
ics by the atom-light coupling. 
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III. MEANFIELD DYNAMICS: THEORY 

The approach we follow in this paper is based upon 
a separation of the field operators into meanfield and 
quantum parts: 

a(t) = a(t) + 5a(t) (5a) 
V(x,t) = VNtp(x,t) + 8V(x,t). (5b) 

In the meanfield approximation the light is assumed to 
be in a classical state with amplitude ait) — (&(£)), 
where |a(i)| 2 corresponds to the average number of pho- 
tons in the cavity, and the atoms are assumed to all 
share the same single-particle wave function ip(x,t) — 
(&(x, t))/\N. The equations of motion for the mean- 
field amplitudes a(t) and ip(x,t) are 



dait) 
i — , — 
dt 

.dip{x,t) 
1 dt 



-A c + NU (cos 2 (x)) - in] a(t) + ir\ (6a) 
' d 2 



dx 2 



\a(t)\ 2 U cos 2 (x)- fx 



tp(x,t) 
(6b) 



where the second equation has the form of a Schrodinger 
equation. The expectation value 



(cos 2 (x))(t) = / dx\ip(x,t)\ 2 cos 2 (x) 



(6c) 



that appears in the first of these equations provides the 
time-dependent coupling between the atomic probability 
density and the cavity mode function. Multiplying this 
integral is the collective atom-cavity coupling parameter 
NUq. When measured in units of the cavity linewidth 
we denote this parameter by /3 



13 = NUo/k . 



(7) 



We illustrate the effect that (3 has on the meanfield dy- 
namics in Figs. [2] and [3] below. 

The time dependence of the meanfield amplitudes due 
to the BOs means that the coupled equations (6a I and 



(6b) already present a nontrivial problem 15 . However, 
in cold atom experiments there is often a large separation 
of time scales between the light and atom dynamics: the 
intrinsic dynamics of the light field is controlled by the 
size of A c and k, whereas the dynamics of the atomic 
meanfield is mainly determined by the period and depth 
s(t) = Uo\a(t)\ 2 of the intracavity lattice. In the experi- 
ments j27J the former is two or three orders of mag- 
nitude faster than the latter. Thus, on the time scales 
relevant to the slow atomic dynamics of the atoms, the 
light field relaxes to an instantaneous "steady state" 



a(t) 



IT] 



NU (cos 2 (x))(t) +ik 



(8) 



which simply tracks the instantaneous value of the atomic 
wave function tp(x,t). This adiabatic elimination of the 



dynamics of the light field results in a clear physical pic- 
ture of the "classical" backaction referred to in the In- 
troduction: As the atomic meanfield evolves in time the 
cavity light field is modulated due to the (cos 2 (x)) term 
in the denominator of Eq. Q which in turn changes the 
amplitude of the light field in the cavity and hence the 
depth of the lattice potential that the atoms feel, and 
so on. Although the adiabatically eliminated form Q is 
quite accurate in the regimes considered here, in our nu- 
merical calculations we nevertheless solve the full coupled 
equations (6a I and (6b). 



In reference [32] we studied the influence the classi- 
cal backaction nonlinearity has upon the band struc- 
ture of atom-cavity systems. The band structure is 
given by the steady state solutions [d = 0, (p(x, t) — 
t pjx ) expj —ifi t/K)] of the coupled equations of motion 
(6a) and (6b) in the absence of the external force /. It 
is straightforward to see that, despite the nonlinearity, 
exact solutions of the steady state problem are given by 
Mathieu functions (like in the linear problem a quantum 
particle in a fixed cosine potential). Mathieu functions 
are Bloch waves and so can be labelled by a band index 
b and quasimomentum q [47] 



fq,b( x ) 



(9) 



where U q j,(x + n) = U q j,{x) has the same period as the 
lattice. In the reduced zone picture q is restricted to lie 
in the first Brillouin zone — 1 < q < 1. Substituting the 
Bloch wave solution into the equations of motion yields 
the steady state equations 



IT] 



NU (cos 2 (x)) + in 



.d_ 

dx 



+ |a ss | 2 [/ocos 2 (a:) 



(10a) 



Uq,b{x) 

(10b) 



where the subscript ss denotes "steady state". Solving 
these equations one obtains a band structure analogous 
to that in the linear case but with the striking differ- 
ence that the nonlinearity can lead to swallowtail loops 
in the bands. It is important to appreciate that this band 
structure is not for the atoms alone, but for the combined 
atom-cavity system. For example, the eigenvalue \i is ac- 
tually a chemical potential rather than the band energy 
(for the underlying energy functional with the light adia- 
batically eliminated see [32 ) , and another difference from 
the linear case is that the lattice depth s = J7o|a ss | 2 is 
not fixed, but instead depends on the values of {6, q}. So, 
for example, the lattice depth changes during a BO as q 
is swept along the band. 

The external force / breaks the spatial periodicity 
and means that Bloch waves are replaced by Wannier- 
Stark states as the steady state solutions of the equa- 
tions of motion (in fact, in any finite system the Wannier- 
Stark states are resonances rather than true eigenstates 
[48) ) . The spatial periodicity can be restored by applying 
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the unitary transformation <p(x, t) = exp(—iftx)ip(x,t) 
which removes the fx term appearing in the hamiltonian 
in the Schrodinger equation ([6bl and introduces a shift 



ft into the momentum operator 



H = 



dx 



2 + s(t) cos 2 (x) - fx 



d 



H = [-i-^+ ft) + s{t) cos 2 (x) 



(11) 



We denote the frame resulting from this transformation 
as the transformed frame (TF), and the original frame as 
the lab frame (LF). 

Let us now consider the dynamics under the influ- 
ence of the force term. We take the initial atomic state 
(p(x,t = 0) = if(x,t = 0) to be a Bloch state in the 
ground band with quasimomentum q — qg. In the adia- 
batic approximation [49 



<p(x,t) w U q (t),b=o(x)exp 



■W) 



/ dt' fi q{t 
Jo 



>),b=0 



(12) 



i.e. the atoms remain in the ground band but the force 
causes the quasimomentum to sweep periodically through 
the first Brillouin zone in accordance with the Bloch ac- 
celeration theorem 



q(t) = Qo + ft 



as can be seen by comparing Eqns. (10b) and (11) 
fact, a careful analysis [50] shows that Eq. (13) holds 



(13) 



In 



even when adiabaticity is broken and interband transi- 
tions are allowed providing these transitions are "verti- 
cal" , i.e. they conserve q. 

This standard approach to BOs remains valid even 
when the lattice depth is modulated in time, as takes 
place in cavities, because amplitude modulation does not 
break the spatial periodicity of the potential and so can- 
not change q [15]. We therefore find that at any later 
time t, the exact atomic meanfield can be expressed as 



<p(x, t) = exp [i(qo + ft)x] U (t) 



(14) 



In general U(t) is in a superposition of bands and so is no 
longer the steady state solution of Eqns. ( |10a[ ) and ( 10b ), 
although it does retain its Bloch form. The advantage of 
the TF is that the quasimomentum is frozen at its initial 
value and we have 



ip(x, t) = exp [iqox] U (t) 



(15) 



so that it is only the spatially periodic function U(t) that 
evolves in time. From the point of view of numerical 
computation this allows us to work with a basis of pe- 
riodic functions (we normalize our wave functions over 
one period of the lattice). At any given time a relatively 
small number of basis functions can accurately describe 
the atomic meanfield state and this greatly reduces the 
numerical effort in the calculation of BOs. Furthermore, 
since the quantum part of the field operators are also 
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FIG. 2. (Color online) Intracavity optical lattice depth 
s(t) = Uo\a(t)\ 2 in units of the atomic recoil frequency ljr 
plotted as a function of time. The curves, which are each 
for a different value of the collective atom-cavity coupling 
parameter NUq/k, were ob tain ed by sol ving the meanfield 
equations of motion Eqns. (6a I and (6b I and illustrate the 



fact that the change in lattice depth over one Bloch oscilla- 
tion increases with NUo/k- In order to maintain a minimum 
lattice depth of 3-Er as NUq/k was increased by changing 
Uq = {1, 3, 5}uo, where uo = 7 x 10~ 3 ljb., we also changed the 
pumping strength as r\ — {30.7, 24.2, 24.3}«;, giving mean pho- 
ton numbers {458,172,117}, respectively. Other parameters 
for this plot are A c = —0.75 n, k = 345 wr, and N = 5 x 10 4 . 
For all the plots in this paper the force is such that the Bloch 
frequency has the value u-g. = wr/4. 



expanded in the same basis as the meanfield, the frame 
transformation also simplifies the numerical calculation 
of the quantum fluctuations as we detail in Section |V| 

By working in terms of Bloch waves rather than Wan- 
nier states, our approach is predisposed towards treating 
wave functions which are localized in momentum space 
rather than coordinate space. This choice is sensible be- 
cause momentum space is a natural setting for BOs as 
is evident from Eq. ( |13[ ). This is also in line with ex- 
isting experiments demonstrating cold atom BOs in free 
space optical lattices [5] [BHS1 [HI US]) where the initial 
state is generally a fairly narrow wavepacket in momen- 
tum space. In this paper we shall therefore restrict our- 
selves to states that are completely localised in quasi- 
momentum (6- function wave packet). Note that the in- 
clusion of direct atom-atom interactions or an external 
trap will break the periodicity in Eq. (11) and induce 



"non- vertical scattering", i.e. scattering between differ- 
ent quasimomentum states. Nevertheless, as explained 
in [41j . in current experiments, e.g. [28], the periodic- 
ity breaking terms are quantitatively smaller and the es- 
sential features can be brought out by the approach we 
adopt. 
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FIG. 3. The lattice depth s(t) in units of the atomic recoil frequency ojr is shown in (a) and its Fourier transform s(uj) is given 
in (b). We have increased the atom-cavity coupling from Fig. [2]to NUo/k = 7.75. At this larger value some fast fluctuations 
on top of the slow BO become visible. Their frequency is dominated by a harmonic at 10 Ub as can be seen in the inset. 



IV. MEANFIELD DYNAMICS: RESULTS 

We now present our numerical results for the meanfield 
dynamics. The initial state at time t = is taken to have 
quasimomentum q = 0, and be given by the solutions a ss 
and£ /o,o(x) of th e meanfield steady state equations [Eqns. 
(10a I and ( |10b )] for atoms in the ground band. This 
state is propagated in time using the meanfield equations 
of motion [Eq. ( [6a] ) and Eq. (6b)]. Our choices for the 
parameter values {Uo,N,t],A c ,k} will be explained at 
the end of this Section. 

Under the action of the external force the atoms be- 
gin performing BOs, which for atoms in extended Bloch 
states gives rise to a breathing motion of the atomic den- 
sity distribution on each lattice site |15j . However, as ex- 
plained in the Introduction, the atomic dynamics are not 
observed directly in the scheme considered here. Rather, 
the classical backaction [which is mediated by the atom- 
cavity coupling integral given in Eq. ( |6c| ] imprints an 
oscillation on the amplitude and phase of the light field 
at the Bloch frequency wb as shown in Fig. [2] (what we 
actually plot in Fig. [2] is the the intracavity lattice depth 
s(t) = U~o\a(t)\ 2 seen by atoms, which is proportional to 
the number of cavity photons \u\ 2 .) The experimental 
signature of the BOs is the photon current transmitted 
by the cavity, and this is given in the meanfield approxi- 
mation by k verta(t)\ 2 , and hence is directly proportional 
to s(t). 

The size of the backaction is controlled by the collective 
coupling (3 — NUq/k, as is apparent from the different 



curves in Fig. [2] As (3 is increased the change in the 
lattice depth over a Bloch period increases and hence the 
visibility or contrast of the BOs as measured by a photon 
detector outside the cavity increases also. We define the 
contrast e as 



0- 



(16) 



Note that in the results shown in Fig. [2j each curve also 
has a different pumping strength r\ in order to maintain 
the same minimum lattice depth of 3£?r . If the lattice 
becomes too shallow interband transition rates (e.g. due 
to Landau-Zener tunnelling around the band edges) be- 
come so high that the atoms effectively fall out of the 
lattice. On the other hand, if the lattice becomes too 



deep the contrast decreases (see Fig. 9a below and also 
Fig. 5 in [22]). A depth of 3E R gives a reasonable com- 
promise. Therefore, although in the rest of this paper we 
will examining the effects of changing the various system 
parameters, we will always maintain the minimum lattice 
depth at 3-Er (the one exception being in Fig. 9a). This 



also allows us to make comparisons between the effects of 
different parameter values upon, e.g. the quantum fluc- 
tuations, whilst keeping the atomic meanfield dynamics 
as similar as possible. 

As the coupling (3 is increased other effects appear 
apart from an increase in the contrast. These effects are 
visible in Fig. [3] (see also Fig. 2 in [15]). In Fig 



3 a 



we 



see that small-amplitude fast oscillations of the lattice 
depth appear on top of the basic BO. These fast oscilla- 
tions have a frequency corresponding to the gaps between 
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bands and can therefore be associated with non-adiabatic 
effects, i.e. the excitation of higher bands [15]. This in- 
terpretation is consistent with the increase in amplitude 
of the fast oscillations with time that can be discerned 
in Fig. 3a Referring to the Fourier transform of s(t) 



plotted in Fig. |3b[ we see that the basic BO dynamics is 
governed by the fundamental ujb and its low lying har- 
monics, whereas the fast oscillations are clustered around 
the tenth harmonic (see inset) and include a continuum 
of frequencies with some peaks at half harmonics. In this 
context it is important to bear in mind that the band 
gaps change continuously in time as q is swept through 
the Brillouin zone and so a range of frequencies is to be 
expected. 

The precision to which wb can be measured in the 
scheme proposed in this paper depends upon the contrast 
[T5] . From the results shown in Fig. [2] it may therefore 
seem that in order to make the most sensitive measure- 
ment possible one should choose f3 to be as large as pos- 
sible. However, this is false for two reasons. One is the 
effect of quantum fluctuations due to measurement back- 
action which is also controlled by f3 and will be the focus 
of Section VIII Another reason, which enters even at the 
meanfield level, is the possibility of bistability in cavity 
photon number for large values of j3 (when the pumping 
is sufficiently large). In [23] this bistability was studied 
experimentally in a uniform unaccelerated condensate, 
which in our language has a quasimomentum q = 0. In 
[3"2"] we studied this problem theoretically and generalized 
it to include finite q: we showed that bistability arises 
from the appearance of swallowtail loops in the bands. 
In the semiclassical picture of a BO the quasimomentum 
scans adiabatically through the entire band and so when 
it encounters a swallowtail loop the system can follow 
a branch that suddenly terminates at some later time, 
leading to fundamentally nonadiabatic behavior |15l 15 1] . 
Hence, in a scheme to measure BOs, it would be bet- 
ter to be in a parameter regime where the cavity is not 
bistable for any value of q. In Fig. [4] we plot the pump 
strength required to maintain the lattice depth at a mini- 
mum value of 3-Er. as a function of (3. The red (solid) and 
blue (dot-dashed) lines enclose the values of rj for which 
the steady state photon number in the cavity displays 
bistability for at least some values of the quasimomen- 
tum. We see that for /3 values as large as 25 (at the fixed 
detuning A c = —0.75k) one can avoid bistability and get 
large contrast in the lattice depth evolution. 

Having emphasized that our choice for the pumping 
strength r\ is guided by the tradeoff between contrast 
and bistability according to Fig. [4j let us now explain 
how we chose the rest of the system parameters used in 
the calculations. There arc three parameters we hold 
constant throughout this paper; the first is the cavity 
damping rate k — 345 which is the value realized in 
the experiment [28] . As a guide to the magnitude of the 
atomic recoil frequency «r used as the frequency unit, we 
note that for 87 Rb atoms in 780nm light ojr — 2tt x 3.8 
kHz. The second constant parameter is the Bloch fre- 
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FIG. 4. (Color online) Plot of pump strength (dashed black 
line) required to maintain a minimum lattice depth of 3£"r 
as a function of NUq/k. The red (solid) and blue (dash dot- 
ted) lines enclose the values of rj for which the steady state 
photon number in the cavity is bistable for any value of the 
quasimomentum of the atomic wave function. One sees that 
for NUq/k ~ 25, the pump strength required to maintain the 
lattice depth leads to bistability. Other parameters for the 
plot are A c = -0.75 k, k = 345 u R , N = 5 x 10 4 



quency wb = Wr/4. The gravitational force on 87 Rb 
atoms in a 780nm lattice provides a Bloch frequency 
very close to this value. Finally, unless specified oth- 
erwise, we keep the meanfield atom number fixed at 
N = 5 x 10 4 and vary Uq in order to vary (3. This last 
choice is motivated by a scaling symmetry of the mean- 
~-d6bb] 



field equations [Eqns. (6a|-(|6bJ], whic h als o hold s for the 
quantum operator equations [Eqns. (24a|-(24b)] below; 
solving the coupled equations for the set of parameters 
{Uo, N, rj, A c , k} is exactly the same as solving them for 
{Uqt, N/r,rj/y/r, A c , k}, where r is some positive scal- 
ing factor. In both the scaled and unsealed versions the 
lattice depth s(t) = t/ |a(i)| 2 is maintained at the same 
value. Thus, each specific calculation one performs repre- 
sents a family of parameters. Our choice for N keeps the 
atomic density dilute enough in a typically sized cavity 
that the approximation of ignoring collisional atom-atom 
interactions remains valid. There is some latitude in the 
choice of A c , but the contrast one obtains at a given 
value of j3 is larger for A c closer to the cavity resonance. 
On the other hand, one also has to make sure that the 
effective cavity detuning 



a: 



a _ 



A c - iVry (cos 2 (:c)) 



(17) 



is less than zero so that we are in the cavity cooling 
regime for the fluctuations [JD] (see Section VI). We set 
A c = — 0.75k since we find that it maximizes the contrast 



for the coupling value of NUq/k — 1. We will examine 
the effect of changing the number of atoms N and the 
minimum lattice depth when we examine the signal-to- 
noise ratio in Section [VI III 



V. QUANTUM DYNAMICS: THEORY 

The approach we take to quantum dynamics is based 
upon a linearization about the meanfield solution, re- 
taining the quantum operators 5a and Sty only to first 
order in the equations of motion. This corresponds to 
the Bogoliubov level of approximation [521 [S3] , suitably 
generalized to describe coupled atomic and light fields. 
A new feature of our problem in comparison to previous 
linearization-based treatments of cavity-QED systems, 
e.g. [33 EU H2] , is that our meanfield is time-dependent 
because of the BOs. This means that the fluctuation 
modes, which must be orthogonal to the meanfield mode, 
also evolve in time (not just their occupations). 

Linearizing about the meanfield solution may appear 
to be an innocent strategy, but, as is well known from the 
theory of Bose-Einstein condensation, care must be taken 
with such U(l) symmetry breaking approaches because 
they introduce a macroscopic (meanfield) wave function 
with a particular global phase at the cost of particle num- 
ber conservation |54) . In particular, when performing a 
linearization about the condensate there is always a triv- 
ial fluctuation mode parallel to it with zero frequency 
(the "zero mode" ) which corresponds to unphysical fluc- 
tuations of the global phase. These issues are even more 
acute when the condensate is time-dependent and the 
boundary between condensate and fluctuation is further 
blurred [55] . 

The zero mode problem can be handled by only in- 
cluding fluctuations that are at all times orthogonal to 
the meanfield. We achieve this by applying the projector 

P(*)|55II551 



so that 

5tyj_(x,t) = P5ty(x,t) (19) 
= / dy[8(x-y)-<p(x,t)<p*(y,t)]6ty(y,t) . 



One consequence of this is that the commutator between 
atomic fluctuations is given by |57j 



5ty ± (x,t),5ty{(y,t)] = (x\P(t)\y) (20) 
= S(x -y)- <p(x,t)ip*(y,t) . 



Unlike, the usual bosonic commutator for the fluctuation 
field 5ty, this is time dependent. 

Next, we transform the atomic fluctuation operator 
from the LF to the TF 

6*{x,t) = 5i>{x,t)e- lftx . (21) 

This simplifies the calculation for the same reasons as 
mentioned in Section |III| for the meanfield. Since only 
"vertical" fluctuations between bands can occur, both 
the fluctuations and the meanfield have the same quasi- 
momentum, which in the TF is frozen at its initial value. 
This allows us to expand the fluctuations and meanfield 
in the same basis 



tp(x,t) = y^c n e 

n 



/2nx 



(22) 
(23) 



where we have set the initial quasimomentum qo = 
without loss of generality. Meanwhile, back in the LF, 
the quasimomentum evolves according to the Bloch ac- 
celeration theorem given by Eq. (13). 



P(t)=I-\cp(t)}(<p(t)\ 



(18) 



We can now write down the coupled equations of mo- 
tion for the fluctuation operators as 



d f ~ 

i — Sa(t) = A(t)6a(t) + VNU a(t) / dxcos 2 (x) <p*(x, t)Sty±(x, t) + <p(x, t)5ty\ (x, t) -Mv^&i) (24a) 

at J i J 

i— 5ty±{x,t) = H(t)Sty±(x,t) + VNU P(t) cos 2 (x)(f(x,t) [a*(t)6a(t) + a(t)6^ (t)] (24b) 



where A(t) = (-A c + NU (cos 2 (x))(t) - in). The struc- 
ture of these equations is such that if initially 8d — 
Sty± = then in the absence of the Langevin terms 
these quantum parts of the fields would remain zero for 
all time. The Langevin fluctuations appear as an inho- 
mogeneous term in the cavity field equation and act as a 
source that drives the growth of da which in turn drives 



the growth of Sty± via the atom-cavity coupling. 

As pointed out in |41j . the dynamics of the complex 
valued operators in the above equations can be solved 
either by separating out their real and imaginary parts 
(optomechanics approach) or by simultaneously solving 
the equations for the hermitian conjugates of the oper- 
ators (the Bogoliubov-de Gennes approach). We choose 
the latter. Collecting the fluctuations into the column 
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vector R{t) = (s&, Sa\ 5*±S& ± ) T , and the noise op- we obtain the operator matrix equation 



erators that act as source terms into the column vector Q „ 

Z(t) = s/2k (1 |t,0, 0) T , where T denotes transposition, *dt R = Mi?(<) + lZ ® ^ 

with 



M(t) 



.4 


VNU a*W(x) 
-VNU Q a*W^{x) 





A* 

NU aW(x) 
NUoaW^x) 



NU aV* 
NU a*V* 
P H(t) 




NU aV 
NU a*V 


pt n(t) . 



(25b) 



where we have introduced the operators 

V • g(x) = J dx (p(x 7 t) cos 2 (x)g(x) 
W(x) = Pcos 2 (x)(fi(x,t) 



(25c) 
(25d) 



i.e. V is an integral operator that acts on a function g(x). 
Since they fall on the off-diagonals, the terms involving 
V and W couple the cavity and atom fluctuations. Ob- 
serve, however, that in the linear approximation used 

here the atomic fluctuation operators S^S±(x, t) are not 
directly coupled to the cavity fluctuation operators 8a(t) 
because this would lead to terms which are of second or- 
der. Rather, the coupling between the two sets of quan- 
tum fields is mediated by the meanfields a(t) and (p{x, t). 

The matrix M(t) is non-normal, i.e. it docs not com- 
mute with its Hermitian adjoint and its left and right 
eigenvectors are not the same. However, it does have the 
following symmetry property: a linear transformation T 
that swaps the first and second, and simultaneously, the 
third and fourth rows, produces a matrix which is pro- 
portional to the complex conjugate of the original [JD] 



T.M.T = -M* 



(26) 



This symmetry, which is a general feature of Bogoliubov- 
de Gennes type equations implies that the eigenval- 
ues (and the associated eigenvectors) occur in pairs of 
the form ±w„ + ij n i.e. with the same imaginary parts 
but with real parts of opposite sign. We shall explore the 
spectrum of the fluctuation matrix M further in the Sec- 
tion IVII We also note that when written in matrix form 
the role of the projection operator becomes clear since 
one can immediately see that the vectors (0, 0, <p(x, t), 0) T 
and (0, 0, 0, (p* (x, i)) span the zero eigenvalue subspace 
of the matrix M and the trivial fluctuations live in this 
subspace. 

The time evol ution of the fluctuation operators is given 
by solving Eq. (25a). However, measurable observables 



are given by expectation values and correlation functions 
of these operators rather than by the operators them- 
selves. To this end we consider the covariance matrix 



C(t) associated with the vector R 

C jk (t) = (RjR k ){t). 



(27) 



Particular cases of Cjk(t) include the total number of 
photonic and atomic fluctuations 



6n(t) = (Stf (t)6a(t)) 



SN(t) 



dx(5^\(x,t)Si'_ L (x,t)). 



(28) 
(29) 



The latter correspond to the number of atoms excited out 
of the meanfield component (i.e. the atomic depletion). 
To obtain the time evolution of the covariance matrix, 



consider the formal solution to Eq. ( 25 ) [58] 
R(t) = G(t,0)R(0) + G(t,0) / G _1 (T,0)Z(r)dT (30) 

where G(t) is a matrix satisfying: 

G(t, 0) = -iM(t, 0)G(t, 0) ; G(0, 0) = 1 . (31) 

We drop the dependence of G on the initial time for 
notational convenience in what follows. Inserting this 
formal solution in Eq. (27) we find 



C(t) = G(t)C(0)G T (t) + G(t)£(t)G T (i) 



(32) 



E(t) 



o Jo 



G-\T){Z{T)Z{T'))[G-\r')YdTdr' 



(33) 



Using the property of the Langevin noise terms given in 
Eq. (4c), we can simplify £(t) as 



E(t) = / G- 1 (t)£»[G- 1 (t)] t c/t 
Jo 



D jk = 2k,8»8. 



'jlOfe2. 



(34a) 
(34b) 



Our main numerical task is thus to solve the matrix dif- 
ferential equation given by Eq. ( |3I| ). In addition, the 
matrix elements of M(i) have to be computed from the 
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meanfields {ct(t), (p(x, t)} obtained by solving the coupled 
equations Eq. ( 6a ) and Eq. (6b I. These latter equations 



are simply a set of ordinary differential equations that we 
solve using an adaptive time-step Runge-Kutta scheme. 
We then solve the matrix differential equation for G(i) 
using the same time grid as the meanfield solution. For 
the matrix differential equation, and the associated solu- 
tion for the covariance matrix C(t), we can again use a 
Runge-Kutta algorithm or exponentiate the fluctuation 
matrix M(i) over the (small) time step intervals [59] . 

As a check on the results we can use the fact that the 
elements of the covariance matrix C(t) have to obey the 
commutator relations Eq. (201 for the operators mak- 
ing up R(t). For example, when the atomic operator is 
expanded as in Eq. (22), the expectation value of the 



commutator relation Eq. ( 20 1 gives 



{8c n 5c\ n ) ~ (Sc^Sdn) = 8 nrn - (n\tp(t))(ip(t)\m). (35) 

The left hand side of this equation gives the difference be- 
tween certain entries of the covariance matrix and we can 
calculate its expected value (the right hand side) from the 
meanfield solution. The degree of agreement between the 
two sides provides a measure of the accuracy of the fluc- 
tuation calculation. In general, we find that the accuracy 
can be increased by taking smaller time steps. 

In closing this section, we would like to point out that 
the meanfield solution already includes Landau-Zener 
type tunnelling that causes the coherent excitation of 
higher bands. By contrast, the effect of Langevin fluc- 
tuations £ is to incoherently populate different bands. 
Within the linear approximation used here, the deple- 
tion of the atomic meanfield by quantum excitations is 
not self-consistent, i.e. the meanfield is always normalized 
to N atoms, whatever the number of depleted atoms 8N. 
The linearized equations used here are clearly only valid 
when 8N << N, as expected from a Bogoliubov-type 
approach. 



VI. SPECTRUM OF ELEMENTARY 
EXCITATIONS 

Before presenting the results of the combined meanfield 
and quantum dynamics (see the next section), we shall 
first examine the excitation spectrum of the atom-cavity 
system. The excitation spectrum gives insights into the 
dynamics and will also be of use in explaining resonances 
that strongly affect the signal-to-noise ratio of the Bloch 
frequency measurement, a topic we will discuss in Section 

m 

We first note that there are two distinct types of ex- 
citation, and hence spectra. The coupled atom-cavity 
band structure discussed in Section (TTTJ refers to mean- 
field excitations which are labelled by a band index and 
a quasimomcntum. They involve every atom and photon 
responding identically since, by the nature of the mean- 
field approximation, they are assumed to be described 



by a single wave function tp(x, t) and the coherent am- 
plitude a(t), respectively. On top of these, there are 
also elementary excitations or quasiparticles whose ener- 
gies are the complex eigenvalues ±oj„ + i"/ n of the matrix 
M(£) given in Eq. ( 25b ). A clear description of the differ- 
ence between the meanfield and the quasiparticle spectra 
for a BEC in a (non-cavity) optical lattice can be found 
in |60) and references therein. In the atom-cavity sys- 
tem the quasiparticles correspond to single quanta of the 
combined fields and are thus polaritons. The fact that 
they come in pairs can be interpreted as an analogue of 
particles and antiparticles |11) . 

Whereas the meanfield band structure is always real, 
the quasiparticle energies have an imaginary part which 
comes from the leaking of the cavity field out of the cav- 
ity. If 7„ < 0, we have dynamical stability and |7„| can 
be interpreted as the lifetime of the quasiparticle. This 
damping effect has potentially very important applica- 
tions in cavity-assisted cooling (3^1 [HT] . If, on the other 
hand, 7„ > we have dynamical instability and heating. 

In general, the elementary excitations have a band 
structure all of their own, i.e. the solutions of the 
Bogoliubov-de Gennes equations take the form of Bloch 
waves with a band index and quasimomentum that can 
differ from that of the meanfield solution about which 



we are linearizing. However, as discussed in Section III 



here we only allow excitations that preserve the quasi- 
momentum (vertical transitions), and thus our quasi- 
particles have the same quasimomentum as their parent 
meanfield solution. Some examples are plotted as a func- 
tion of quasimomentum in Fig. [5] We note in passing 
that there is in fact no need to quantize the fluctuations 
5^ and Sa about the meanfield in order to obtain M. 
Thus, the spectrum of elementary excitations can also be 
obtained by considering classical fluctuations about the 
meanfields. This is well known in the case of BECs, where 
the linearized Gross-Pitaevskii equation (the meanfield 
equation of motion) gives the Bogoliubov spectrum. The 
truly quantum part of the dynamics here comes from the 
Langevin fluctuation operators £(i). 

The eigenvectors of M can be roughly classified into 
three kinds: cavity-like modes, hybridised atom-cavity 
modes, and marginally stable modes |40j . The cavity- 
like modes (depicted by the green dash-dotted lines in 
Fig. [5| are close to being pure cavity field modes with 
only a small atomic component. Hence, their eigenval- 
ues have a real part with magnitude close to the effec- 
tive detuning A° = A c — iVL 7 o(cos 2 (cc)), and an imag- 
inary part approximately equal to — n. The hybridized 
modes (depicted by the red solid and blue dashed lines in 
Fig. [5]) have some atomic and some cavity field properties, 
whereas the marginally stable modes are purely atomic in 
nature with zero cavity component. As we shall demon- 
strate below, the marginally stable modes occur at the 
points q = and q — ±1, i.e. at the band center and 
edges, and their name derives from the fact that their 
imaginary part is zero. 

The properties of the hybridised and marginally stable 
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(a) Real part of the quasiparticle spectrum as a function of 
quasimomentum 
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(b) Imaginary part of the quasiparticle spectrum as a function 
of quasimomentum 



FIG. 5. (Color online) Low lying levels in the quasiparticle spectrum (elementary excitations). The frequency of the nth level 
is generally complex uj„ + i"1n- The parameters used in the plots are Uo = 0.01 uj-r, k = 345 ujr, A c — —0.75 k, and N — 5 x 10 4 . 
The red (solid) and blue (dotted) lines correspond to hybridised atom-cavity modes and generally have non-zero imaginary 
parts except at certain special points such as at the band center and edges where they can become marginally stable and 
decouple from the cavity The green (dash dotted) line corresponds to a cavity-like mode, i.e. the real part of its frequency is 
close to the effective detuning frequency A° ff (q), and the imaginary part is close to —k. 



modes are determined by the sign of A® . When A° ff < 
we find 7„ < and we are on the cooling side of the effec- 
tive resonance. On the contrary, when A° ff > we find 
7„ > and we are on the heating side. A calculation 
of the dynamics on the heating side is not stable since 
the linearisation will fail after a short time due to the ex- 
ponentially growing number of quasiparticles. Thus, the 
calculation of the spectra serves a very useful purpose: it 
guides our choice of A c so as to ensure that we are always 
on the cooling side of the resonance. 



In Fig. 5a the red (solid) and blue (dashed) lines give 
the magnitudes of the real parts of the frequencies of 
the two lowest eigenmodes as a function of quasimomen- 
tum. The magnitudes of the imaginary parts are plot- 
ted in Fig. [5b] Notice that the imaginary part of one 
of the modes goes to zero at the band edges (red solid 
line) and the other goes to zero at the band center (blue 
dashed line) . This implies that these excitations are only 
marginally stable at those specific values of the quasi- 
momentum. The vanishing of the imaginary part of the 
frequencies at these points can be understood as follows: 
a Bloch wave with q = (q = 1) is even (odd) about the 
center of a single cell [0, it] of the cos 2 (a;) potential. Since 
the atomic meanfield solution (p(x) is a Bloch wave it has 
well defined parity at these points. The same is also true 
for the eigenmodes 8^/±(x) of the fluctuation matrix M, 
for these are also Bloch waves. In this case the integral 



f dxtp* (x) cos 2 (x)S*S> ±(x) will sometimes vanish identi- 
cally because the integrand can contain functions with 
opposite parity. Examination of the fluctuation matrix 



M given in Eq. (25b) shows that it is exactly this in- 



tegral that controls the weight of the cavity part of the 
quasiparticle eigenmode, and so at q — 0, ±1 we can have 
undamped quasiparticles with 7 = 0. For nonzero values 
of quasimomentum the meanfield wave function has no 
particular parity and there are no marginal modes. 



VII. QUANTUM DYNAMICS: RESULTS 

In this section we present results from the numerical 
solution of the quantum equations of motion. We as- 
sume that at t = the fluctuation fields 8a and 8^± are 
in their vacuum states and expand the atomic part in 
the basis given in Eq. (23). This allows us to construct 
the covariance matrix Eq. (27) C(t = 0) which we then 
evolve to later times using Eq. (31). In order to per- 



form this task we need the fluctuation matrix M(i) as 
a function of time which in turn requires the meanfield 
solution {(p(x, t), a(t)} as input. We must therefore solve 
the meanfield dynamics on the same discretized time grid 
in parallel to the computation of Eq. (31 ). Once we have 



computed C(t), we can extract the various physical quan- 
tities of interest from it. 
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(a) Excited atom fraction versus time for NUq/k = 0.1 (b) Excited atom fraction versus time for NUq/k = 1 

FIG. 6. (Color online) Growth of the excited atom fraction over five BO periods for weak (a) and moderately strong (b) 
atom-cavity coupling. The red (solid) curves are given by solving the full quantum problem in the form of Eq. ( |31[ ), whereas 
the black (dashed) curves are the result of making the approximation that the quantum fluctuations come purely from vacuum 
shot noise, i.e. the coherent state approximation. The meanfield dynamics for (b) is given by the red (solid) curve in Fig. [2] 



One such quantity is the number of atomic excitations 
SN, which is defined in Eq. (29). This should not be 



confused with the number of quasiparticles because the 
latter are generally made up of both atomic and cav- 
ity field components. If it were not for the Langevin 
noise, the evolution would be perfectly coherent and SN 
would be zero. However, the presence of Langevin noise 
in the electromagnetic field generates atomic excitations 
via the atom-cavity coupling. In Fig. [6] we plot the frac- 
tion SN(t)/N as a function of time for five Bloch periods 
for two different coupling values. The red (solid) curves 
are given by a full solution of the quantum equations, 
whereas the black (dashed) curves are made with a co- 
herent state approximation for the cavity field, which will 
be outlined below and is discussed in more detail in Ap- 
pendix [A] The gradient of these curves gives the heating 
rate and over five BOs this appears to be roughly linear. 
However, over longer times the number of atoms (and 
photons) excited out of the meanfield begins to saturate, 
as shown in Fig. [7} The saturation occurs because on the 
cooling side of the effective resonance, where A° ff < 0, 
there is a competition between the heating due to the 
Langevin noise Z and the dynamical stability of M lead- 
ing to nonzero but finite values of SN at long times. 

Superimposed upon the steady growth of the number 
of atomic excitations shown in Fig.[6j it is possible discern 
small modulations with a period equal to that of a BO. 



These modulations are much more evident in the photon 
dynamics shown in Fig. |7b| They are due to the oscil- 
lation of the atom-cavity coupling during a BO, which 
takes its greatest value at q = and its smallest value at 
q — ±1. Indeed, close examination of the modulations 
in Fig. [6] reveals that the heating rate at the zone edges 
is lower than at the zone center (recall that the initial 
condition is q — at t — 0, and the system returns to 
this point after time Tb: at time Tb/2, the system is 
at the edge of the Brillouin zone). The coupling inte- 
grals enter M, and referring to Fig. [5j we that this effect 
manifests itself in the dependence of the quasiparticle 
mode energies upon q: at the zone edges the quasipar- 
ticle mode with the smallest real part (red solid curve) 
becomes marginally stable, i.e. the cavity light field part 
and the atomic part decouple. This means that close to 
the zone edge new quasiparticles can only be created in 
a mode with higher energy (blue dotted line) which leads 
to a lowering of the heating rate. 

In order to gain further insight into the fluctuation dy- 
namics, let us develop a simple approximate model which 
we will call the "coherent state approximation" that as- 
sumes that the cavity field is described by a Glauber co- 
herent state. This is a quantum state that has the same 
fluctuations as the vacuum and so neglects the fluctu- 
ation correlations that develop between the atoms and 
cavity field. Consider the exact solution to the first or- 
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FIG. 7. (Color online) Plots of atomic (a) and photonic (b) fluctuation occupation number over 20 BOs calculated using a 
numerical solution of Eq. ( |31[ |, In (b) the lowest curve (red) is for j3 — 1, the middle curve is for /3 = 3 and the highest curve 
is for P = 5. Note the saturation in the growth rate, which is especially obvious in the photons. The inset in (b) shows a close 
up of the photonic fluctuation number as a function of time for /3 = 1. It can be seen how after a transient period the photonic 
fluctuation number oscillates at the Bloch period, thereby mirroring the meanfield dynamics. 



der inhomogeneous differential equation for cavity field 
fluctuations Eq. ( 24a ) , which can be formally written as 



In our earlier discussion (Section IV I of desirable param- 



Sa(t) 



-i /* dt' A(t') , 



(36) 



dt'e^o dt '' A ^')[V^{t')-iVNU a{t')X{t')] 



where X(t) = J dx ip* (x, t) cos 2 (x)5'&±(t) + h.c, and 
A(t) — — A° ff (i) — in, as above. From Eq. (36) we see 



that the cavity field fluctuation has two distinct contri- 
butions: the first term depends on the Langevin noise 
which accounts for vacuum fluctuations, whilst the sec- 
ond term depends on the state of the atoms. The coher- 
ent state approximation consists of dropping the latter 
term in favour of the former to give 



Sa(t) « / dt'e^^-^-^iit') . 



(37) 



In writing 5d(t) in this way we have taken advantage 
of the fact that the cavity decay rate k is much faster 
than the frequency ub at which A° s (t) evolves, and so 
the integrand is appreciable only for times t — t' < k^ 1 
during which A° is a constant and can be evaluated at 
time t. 

The regime of validity of the coherent state approxima- 
tion can be estimated from its derivation which requires 
r = y/NU \a(t)\/V2K « 1. Note that r 2 = /3s(t)/2. 



eters, we stipulated a minimum lattice depth of s(t) ~ 
3wr, which implies that the validity of the coherent state 
approximation here is contingent upon j3 -C 1, i.e. this is 
a weak coupling approximation. 

A measure of the magnitude of the correlations be- 
tween the atomic and cavity field excitations, i.e. the de- 
gree to which the cavity field differs from a coherent state, 
can be seen in the difference between the red solid and 
black dashed curves in Fig. [| for SN(t) /N. The latter 
neglects these correlations. We can thus infer that when 
j3 = 0.1 the correlations are small, but when (3 = 1 their 
effect is noticeable, although not enough to make a qual- 
itative change over 5 BOs. In fact, as mentioned in the 
Introduction, the heating rate of a cloud of cold atoms 
inside a cavity has been measured in the experiment by 
Murch et al |31j , and they found it to be consistent with 
the predictions of vacuum noise. The new feature in our 
problem is that the effective cavity drive detuning A° ff (t), 
which appears in the phase terms in Eq. (37), is Bloch 
periodic due to the meanfield dynamics. 

The assumption of vacuum noise is often made in cav- 
ity optomechanics and allows some results to be derived 
analytically [3DJ l6"2"r[6"5] . The paradigmatic example of a 
cavity optomechanical system is a cavity with one end 
mirror attached to a spring or cantilever i.e. a harmonic 
oscillator driven by radiation pressure. Although ultra- 
cold atoms in a cavity can sometimes be mapped onto 
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this problem [Ml 133 EH ESI M\ , that is not the case for 
the present system because the atomic Bloch states do 
not map very faithfully onto a single harmonic oscillator. 
Nonetheless, inspired by the optomechanical example we 
can tackle the problem by mapping it onto a collection of 
independent oscillators. The atomic excitation occupa- 
tion 8N(t), which is plotted as the black (dashed) curves 
in Fig. [6j is then the sum over the occupation numbers of 
these independent oscillator modes 5N(t) — SNj(t). 
The details of this procedure will not be described here 
but are presented in the Appendix[S] Suffice it to say that 
the basic idea is to consider the noise as a perturbation 
to the oscillator dynamics and then use Fermi's golden 
rule to develop a rate equation describing the occupation 
number dynamics of each oscillator |63| 



d(5Nj 
d~T 



II.) 



(38) 



which is Eq. ( A8 ) in Appendix [A] In this expression r„j 
and r^- are the transition rates "up" and "down" for 
the jth oscillator and they are proportional to Sj^j^(—u)j) 
and Sjrjr(ujj), respectively, where Sjf„f(w) is the spectral 
density of force fluctuations (shot noise power spectrum) . 
Thus, each oscillator is driven and damped by vacuum 
noise, with the rates of driving and damping being time 
dependent (due to the meanfield BO dynamics). Such a 
description allows an equilibrium to be reached and is the 
origin of the saturation effect visible in Fig. [7] Indeed, 
setting the left hand side of the rate equation to zero we 
obtain an expression for the equilibrium occupation of 
each oscillator 



r dj (t)-r uj (t) 



(39) 



By equilibrium here we mean a dynamic state in which 
SNj is slaved to the instantaneous values of the transition 
rates. In the next two sections we examine the implica- 
tions of the fluctuations for a precision measurement, i.e. 
how the fluctuations put a limit on how large a value of 
/3 can be chosen for a precision measurement. 



VIII. SIGNAL-TO-NOISE RATIO: THEORY 

We now explore how the inclusion of quantum noise 
affects the precision measurement proposal in |15) . Recall 
the basic idea shown schematically in Fig. [I] a cloud 
of cold atoms undergoes Bloch oscillations (e.g. due to 
gravity) inside a Fabry-Perot cavity, and the light field 
transmitted through the cavity is measured in order to 
determine the Bloch frequency. In order to quantify the 
measurement performance we will compute the signal-to- 
noise ratio using standard input-output theory [55]. 

Let us consider a double sided cavity with mirrors with 
matched reflectivities providing equal amplitude damp- 
ing rates of n/2. The quantum part of the input fields for 
both the top (driving side) and bottom (detection side) 



mirrors is given by the electromagnetic vacuum. Since 
we are not going to consider classical fluctuations of the 
driving laser we do not include the meanfield part a of 
the input field, but introduce it via the hamiltonian in 
Eq. ([3| . In our consideration of system dynamics in ear- 
lier sections we implicitly assumed a single sided cavity 
giving an amplitude damping rate of k, and associated 
with this decay is a Langevin noise term \/2K^(t). In a 
double sided cavity we have two independent noise terms 
of the form v^{£* CO j It can be shown that the 

dynamics of the intracavity system (both meanfield and 
fluctuations) are independent of whether we assume a 
double sided or single sided cavity as long as we divide 
the net damping equally among the two mirrors (pro- 
vided they have matched reflectivities) . The transmitted 
light field is the output field at the bottom mirror which 
is related to the input field at the bottom mirror as 

flout (0 = -a- m (t) + y/Ka(t) = -i b (t) + y/iia(t) (40) 

where ai n and a out in this equation refer to the fields 
at the bottom mirror. The transmitted photon current 
is given by the operator I ut(t) = al ut (t)a ou t(t), where 
again a out refers to the field leaving the bottom mirror. 

An experimentally straightforward method for measur- 
ing the Bloch frequency consists of recording the trans- 
mitted photon current using a photodetector. It is useful 
to consider the Fourier transform of the data [69] 

N(cj,T)= [ dtcas(u>t)I ovlt (t) (41) 
Jo 

and define the signal-to-noise ratio for the measurement 

as 



|(^,T))| 2 
SNR= AN^,T) 



(42) 



where AN 2 {w,T) = ((N-(N)) 2 ). Thus, the SNR is the 
ratio of the spectral density of the photon current to its 
variance and provides one measure of the sensitivity of 
the scheme. 

Let us first evaluate the SNR for a classical cavity field 
a(t) = a(t). In this case one finds that the signal ampli- 
tude and variance are given by 

(N(u,T))=k[ dtcos(ojt)\a(t)\ 2 (43) 
Jo 

r T 

AJV 2 (w,T) = k / dt cos 2 (ut)\a(t)\ 2 . (44) 
Jo 

In order to obtain an approximate magnitude for the 
SNR we further assume that the detection rate goes as 
« i?(l+ecos[wB£|) [H]: where e is the contrast parameter 
denned in Eq. (16). Setting the classical photon current 



K|a(t)| in the above formulae equal to this detection rate 



gives 



SNR(w B ,T) 



e 2 RT 



(45) 
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FIG. 8. (Color online) Plots of the SNR as a function of coupling strength j3 (a), and as a function of integration time for 
different values of ft (b). In (a) the SNR was computed for an integration time of 10 Bloch periods (Tb) and the red (dots) curve 
gives the meanfield dynamics plus detector shot noise result, whilst the black (crosses) curve includes measurement backaction, 
i.e. the effect of quantum fluctuations upon the coupled atom-cavity dynamics. In (b) the red (solid) and blue (dotted) curves 
lie almost on top of each other and correspond to values of /3 just before the first dip in the SNR shown in Fig. |8a| whereas 
the green (dash-dotted) curve corresponds to a value of /3 in the dip. For all plots the minimum lattice depth was 3Er,. Other 
parameters are the typical ones mentioned in the text. 



Despite appearances, this result does include quantum 
noise to a certain degree because without the Langevin 
operators the variance given in Eq. p4| would have been 
zero, i.e. even when the cavity field is classical the output 
field contains a quantum part a out = a(t) — Thus, 
the above calculation includes detector shot noise, also 
known as measurement imprecision [53] , but neglects the 
effect of quantum fluctuations on the coupled dynamics 
inside the cavity, i.e. quantum measurement backaction. 
Note that this is a different approximation from the co- 
herent state approximation used in Section |VII[ where 
quantum fluctuations were included in the cavity dynam- 
ics by using a Glauber coherent state, i.e. a state with 
vacuum noise, as cavity field, albeit one whose fluctua- 
tions are unaffected by the presence of the atoms. 



The SNR given by Eq. ( 45 ) predicts that the sensitivity 



of the scheme can be increased indefinitely by increasing 
the mean total number of photons collected RT and also 
the contrast e. The former effect is the standard one ex- 
pected from the general theory of measurements with un- 
corrected fluctuations. The latter is intuitively plausible 
too, but, however, can not be the whole truth because, 
as stated above, it neglects the effect of measurement 
backaction upon the dynamics which is expected to be- 
come important at larger values of j3. When fluctuations 
are included a(t) — a(t) 4- Sa(t), and the mean signal 



amplitude is given by 

r T 

(N(u,T)) =k I dtcos(cot) (\a(t)\ 2 + (5a f (t)8a(t)}) . 
Jo 

(46) 

In fact, this is not so very different from the meanfield 
photon number given by Eq. ( 43 ) because we are by de- 



sign working in a regime where the meanfield dominates 
the fluctuations. However, the same is not true of the 
signal variance. The expression for the signal variance 
including fluctuations is cumbersome and is presented 
in Eq. (Bl| in Appendix fBl For present purposes it is 



enough to note that it includes a collection of terms 
that depend on integrals over two-time correlations of 
the photon fluctuations. These two-time correlations are 
challenging to evaluate numerically not only because the 
fluctuations occur on time scales k _1 much shorter than 
the BOs, but also because they require the storage and 
manipulation of data at two times. Furthermore, the con- 
tinuous driving by the BOs means that the correlations 
are not stationary in time, i.e. they do not just depend 
on t\ — <2, and this forces us to calculate the SNR in 
parallel to the system dynamics starting at t = 0. Unfor- 
tunately, due to limited computing power, we have only 
been able to track the SNR over ten Bloch periods which 
is certainly shorter than the coherence time of the BOs 
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(a) Signal-to-noise ratio as a function of min. lattice depth (b) Signal-to-noise ratio as a function of atom number 



FIG. 9. (Color online) Plots of the SNR as a function of minimum lattice depth s(t = 0) (a), and as a function of atom 
number N (b). In both plots the red (dots) curves were computed from meanfield theory plus detector shot noise, and the 
black (crosses) curves were computed including quantum measurement backaction. For all points NUo/k — 1 and the signal is 
integrated over 10 Bloch periods. In (a) it is evident that the SNR decreases in both cases for larger lattice depths. In (b) it 
is evident that the SNR increases linearly as a function of N in both cases. 



for the parameters we use. An actual experiment would, 
of course, not suffer from this limitation and would bene- 
fit from running until the BO coherence time is reached. 
The main steps of our algorithm for calculating the two 
time correlations are provided in Appendix [Bj 



IX. SIGNAL-TO-NOISE RATIO: RESULTS 

We now show how the SNR depends on the various 
system parameters. Due to the size of parameter space, 
this will not be an exhaustive study, but rather an ad hoc 
choice that nevertheless we hope is experimentally rele- 
vant. We begin by looking at the SNR as a function of the 
coupling parameter — NUq/k. In Fig. 8a we plot the 



SNR evaluated at ujb for an integration time of 10 Bloch 
periods. We change j3 by increasing Uq but also change 
•q to maintain the same minimum lattice depth of 3En 
throughout. The results without measurement backac- 
tion (i.e. the dynamics in the cavity is purely meanfield) 
are plotted by the red (dots) curve which monotonically 
increases until about /? = 12. The initial increase of the 



SNR with (3 is in line with expectations based on Eq. ( 45 1 



The turnover of the red curve near f3 = 12 is in a sense 
an artefact that arises from having evaluated our SNR at 
ujb: it so turns out from the meanfield solution that for 
j3 > 12 the fraction of the power in the fundamental of 
s(uj) begins to decline and is diverted to higher harmon- 
ics. However, there is no real reason other than simplicity 



to only consider SNR(wb) (any harmonic of wb gives in- 
formation about the applied force and inclusion of all of 
them in the data analysis would extract the maximum 
possible information from the measurement). The full 
calculation including measurement backaction is plotted 
by the black (crosses) curve. The first thing to notice 
is that measurement backaction always lowers the SNR. 
Secondly, the full SNR monotonically increases only until 
ft « 7, and thereafter suffers from dramatic dips which 
we explain below as being due to resonances with quasi- 
particle excitation energies. These two observations are 
the main results of this paper. In Fig. [8b] we plot the 
SNR as function of the total integration time for three 
values of j3, two before the first dip in the SNR and one 
in it. This plot further illustrates that for f3 > 7 there is 
a dramatic lowering of the SNR. The BO dynamics are 
also clearly visible due to the fact that the contrast is 
periodically growing and shrinking as the lattice depth 
grows and shrinks. 

integration time of 10Tb ■ Whereas for the value (green 
solid line) to the right of the peak one can see that the 
SNR has more or less fiat lined near t as 10Tb. Thus 
in Fig. [8a] (with integration time of 10Tb), for points to 
the left of the peak, the SNR is still gaining due to the 
increased measurement strength whereas this is not the 
case to the right of the peak. 

In Figs. [9a] and [9b] we show how the SNR depends 
on other parameters, namely the lattice depth and total 
number of atoms. In particular, in Fig. 9a we plot the 
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quasimomentum (units of k c ) NUq/k 

(a) Lowest excitation frequency for three different values of /3 (b) Lowest excitation frequency for all q as a function of f) 

FIG. 10. (Color online) Plots of the lowest quasiparticle excitation frequency uui about the adiabatic solution. In (a) this is 
given as a function of quasimomentum for three different values of j3: For small /3 (red dashed curve) the minimum of the 
frequency occurs at q — ±1, but for larger values of /3 the minimum shifts in to smaller values of q. Since each quasiparticle 
excitation has an energy varying with q, in (b) we plot the range of possible excitation frequencies contained in cji for all q as 
a function of /3. The frequency units are the Bloch frequency uis. 



SNR as a function of the minimum lattice depth in the 
cavity for the coupling value NUq/k = 1. The lattice 
depth is changed by increasing r\. The red (dots) curve 
gives the SNR calculated using only meanfield dynamics 
plus the effect of shot noise at the detector and justi- 
fies the claim made in Section [TV] that for larger lattice 
depth the contrast decreases. The SNR calculation in- 
cluding measurement backaction fluctuations is given by 
the black (crosses) curve, and has the same qualitative 
behaviour but is somewhat lower. Fig. [9b]plots the SNR 
as a function of N, where for different values of N we 
keep NUq/k = 1 constant by scaling Uq. We also scale 
the pump strength rj to maintain the same intracavity 
lattice depth s(t) in all the cases. As we pointed out in 
Section |IV[ this method of scaling the system variables 
leaves the form of the meanfield and fluctuation equa- 
tions unchanged. The only quantitative change is that 
the meanfield cavity field solution a(t) is scaled by the 
same yfr factor as the pumping. This leads to a linear 
scaling of the SNR as a function of N (with and without 
fluctuations) as shown in the plot. It is interesting to note 
that the rate of increase is different for the calculation in- 
cluding fluctuations compared to that without. Clearly 
there is a gain in the SNR with N, but one has to also 
keep in mind that in our treatment we have completely 
ignored the atom-atom collisional interaction which will 
eventually become important at high atom densities. 

Finally, we shall explain the physical origin of the com- 
plicated series of dips in the SNR when j3 > 7 that are 



seen in Fig. 8a Consider the spectrum of quasiparticle 



excitations about the adiabatic meanfield solution intro- 
duced in Section |VI| For the example shown in Fig. [5a] 
the smallest excitation frequency occurs at the band edge 
q = ±1 and the largest at q = 0. As f) is increased in 
the usual manner (holding the minimum lattice depth 
constant), the g-dependence of the quasiparticle spec- 
trum e volve s, as shown for the quasiparticle mode lu\ 
in Fig. 10a Thus, the range of frequencies (i.e. across 



the entire Brillouin zone) contained in uj\ also evolves 
with f3 and is shown in Fig. |10b| If the meanfield dy- 
namics happens to contain any frequencies that fall in 
this range there is clearly the possibility of a resonance, 
exciting quasiparticles and lowering the SNR. This is ex- 
actly what happens as can be seen from Fig. [IT] which 
plots the total power in the harmonics of wb that fall in 
the frequency range covered by uj\ . The two peaks in Fig. 
[TTIat f3 sa 8 and j3 « 12 coincide exactly with the dips in 



Fig.[8a[ Referring back to the inset in Fig. 3b which was 
deliberately evaluated at /3 = 7.75 for this very purpose, 
we can see the part of the meanfield spectrum that falls 
in the range spanned by u%. In the absence of BOs the 
quasiparticle excitation uj\ is very narrow, with a width 
given by the imaginary part 71 evaluated at q = 0. How- 
ever, the BO dynamics effectively broadens the resonance 
by orders of magnitude to that shown in Fig. 1 10b] and this 
has a dramatic effect on the SNR. 
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FIG. 11. The normalized power in the harmonics of ujb as 
calculated from the Fourier transform of the meanfield solu- 



tion (see Fig. 3b I that lies in the f reque ncy range of the lowest 
quasiparticle excitation (see Fig. 10b I. 



X. SUMMARY AND CONCLUSIONS 

In this paper we have extended our previous analysis of 
Bloch oscillations of ultracold atoms inside a cavity to in- 
clude the effects of quantum noise in the electromagnetic 
field. The quantum noise originates from the open nature 
of the cavity and can be interpreted as a form of quantum 
measurement backaction because it perturbs the dynam- 
ics. The magnitude of the backaction is controlled by the 
dimensionless atom-light coupling parameter (3 — NUq/k 
and we find that it can strongly affect the sensitivity of a 
measurement of the Bloch oscillation frequency wb and 
therefore the determination of the magnitude of the ex- 
ternal force F driving them. 

Our treatment is based upon the coupled Heisenberg 
equations of motion for the atoms and light which we 
linearize about their meanfield solutions, i.e. a Bogoli- 
ubov level approximation. We solve the meanfield level 
dynamics exactly and hence coherent effects such as 
Landau-Zcncr tunneling between bands are fully taken 
into account. A spectral decomposition of the meanfield 
solution shows that it is dominated by wb and its first few 
harmonics, but as (3 is increased spectral power begins to 
spread to higher frequencies. 

Quantum noise is introduced via Langevin operators 
which act as inhomogeneous source terms in the Heisen- 
berg equations. These terms lead to both the creation 
and annihilation of quasiparticles (quantized excitations 
with a mixed atom-light character) and hence can both 
heat and cool the system. If the system is started off with 
no quasiparticles, i.e. at zero temperature, the number of 



excitations initially grows in time at a rate that increases 
with p. However, the number of excitations eventually 
saturates due to competition between the cooling and 
heating processes (provided we are in the cooling regime 
Af = A c - NU (cos 2 {x)) < 0) with the "equilibrium" 
number of excitations being dynamic in the sense that it 
is modulated by the Bloch oscillation. 

In order to gain some insight into the numerical cal- 
culations we used Fermi's golden rule to develop a semi- 
analytic model for the heating rate in terms of a simple 
rate equation for the number of atomic excitations. In 
so doing we approximated the cavity light field by a co- 
herent state whose quantum fluctuations are the same 
as those of the vacuum. This a common approximation 
in cavity optomechanics but ignores the quantum corre- 
lations that build up between the atoms and the light. 
Comparing this with the exact numerical results for the 
number of atomic excitations, we infer that the field is 
close to a coherent state for small f3, but differs from it 
as (3 is increased, as expected. 

The above calculations can be applied to the estima- 
tion of the signal-to-noise ratio for a continuous mea- 
surement of wb- For example, we find that the SNR de- 
creases with intracavity lattice depth, and increases with 
the number of atoms. Our principal result, however, con- 
cerns the dependence upon /3. We find that the SNR can 
be severely reduced due to resonances between the quasi- 
particle spectrum and the Bloch oscillating meanfield for 
certain ranges of f3. Indeed, the SNR behaviour depicted 
in Fig. [8] is much more complicated than that found in 
the standard example of a quantum limited position mea- 
surement of a harmonic oscillator, e.g. the end mirror of 
a resonant cavity [64]. In that system, the SNR is de- 
termined by the competition between the "measurement 
imprecision" (detector shot noise), which decreases with 
increasing measurement strength, and the measurement 
backaction, which increases with increasing measurement 
strength, and correlations between the two can be ig- 
nored to a good approximation. This leads to a smooth 
curve (see Fig. 5 on p. 1171 of the review [64]) with a sin- 
gle maximum at the measurement strength where the two 
effects are equal. This is where the measurement should 
be performed for maximum sensitivity. By contrast, in 
our case we have a cloud of atoms occupying Bloch states 
in an optical lattice and thus our system does not corre- 
spond very well to a single harmonic oscillator (except in 
the limit where the lattice is extremely weak so that the 
atoms are predominantly in a state which is uniform in 
space |28j . but then Landau-Zener tunneling will be so 
severe that the atoms will quickly fall out of the lattice 
when the external force F is applied). Add to this the 
fact that our system is driven by an external force and so 
scans through the entire Bloch band in a time-dependent 
fashion, leading to the possibility of resonances, and it 



is not surprising that our resulting SNR in Fig. 8a does 
not have a simple maximum as a function of (3. However, 
we can make the parameter dependent statement that it 
seems safest to choose (3 < 7 which lies below the point 
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where the resonances set in (and for ft > 25 we find op- 
tical bistability which will destroy the Bloch oscillations 
|32|). The resonances only occur in the calculation when 
quantum measurement backaction is included and so pro- 
vide a salutary example of when the latter is important. 
Nevertheless, away from the resonances the SNR for this 
continuous measurement is large and is in pretty good 
agreement with an approximate calculation based upon 
purely meanfield dynamics in the cavity plus detector 
shot noise. 

In comparison to previously studied cold atom cavity- 
QED systems, or even cavity optomechanical systems, a 
new feature of our Bloch oscillating system is the time- 
dependence of the meanfield. Apart from the resonances 
discussed above, this also has implications for the com- 
putational scheme we use to calculate the results. For 
example, all the fluctuation modes should be orthogo- 
nal to the meanfield mode as well as to each other, and 
hence they also evolve with time. Furthermore, the two- 
time correlation functions that are needed to calculate 
the signal variance that enters the SNR are not station- 
ary in time, meaning that a large amount of data must 
be stored. This is especially true because the Bloch pe- 
riod is roughly three orders of magnitude larger than the 
quantum fluctuation timescale 1/k and hence the cal- 
culation of the SNR over even a few Bloch periods is 
quite intensive in the regime where the coherent state 
approximation breaks down. In non-cavity BO exper- 
iments it has been shown that coherent dynamics can 
run for thousands of Bloch periods [9]. In a continuous 
measurement scheme, such as that proposed here, the 
quantum measurement backaction reduces the coherence 
time but unfortunately we have been unable to go much 
beyond ten Bloch periods with our numerical computa- 
tions of the SNR and thereby find this coherence time for 
our scheme (we have, however, given an estimate in [15] 
based upon the idea that the spontaneous emission rate 



sets the upper limit on coherent dynamics). Nonethe- 
less, our short time calculations illustrate quantitatively 
that it may be advantageous to remain at small j3 and 
integrate for longer times. In future work we hope to de- 
velop approximations and/or better numerical methods 
to treat large values of the coupling. 

frequency and the minimum of the atomic fluctuation 
spectrum leads to enhanced occupation of the fluctuation 
modes. 
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Appendix A: Coherent State Approximation 

In this appendix we provide details of how the coher- 



ent state approximation described in VII can be used to 



derive a simple approximation for the atomic fluctuation 
dynamics. 



Consider Eq. p4bl| for the atomic fluctua- 
tion dynamics, and let us decompose the atomic fluctua- 
tion operator (working in the TF) into the instantaneous 
eigenbasis \v k (t)) of the time-dependent meanfield hamil- 
tonian H(t) 

*#x(t) = J>*(*)>*S fc (t) (Al) 
k 

H(t)\v k {t)) = E k {t)\v k {t)) (A2) 

where we have dropped the coordinate dependence of the 
atomic fluctuation operator and we use bra-ket notation 
for the single particle states. On substituting the decom- 
position Eq. ( Al | into Eq. ( 24b ) we obtain 



dSbj{t) 
dt 



(A3) 



iy/NU (a*(t)6a(t) + a(t)6tf(t)) (vj(t)\P(t) cos 2 {x)\<p(t)) . 



Clearly in the above form the dynamics of the 5b j (t) are 
coupled amongst themselves. To decouple the dynamics 
consider the second term (with the time derivative) on 
the right hand side of Eq. (A3). It can be shown that 
(EDI): 



(vJt)\-\v k (t)) ] t k 1 — 

\ j\ )\ dt \ k\ ii E k (t)-E;j(t) 



to(t)I^M*)) 



Now for k = j ± 1, the matrix element goes to zero for 
symmetry reasons and for other k, the energy difference 



will make the term smaller. Thus we need to only take 
care of the term for k — j. Focusing on that now we have 



h.c. = 



i.e. the derivative term with j — k is purely imaginary. 
For the time dependent hamiltonian H(t) the potential 
term cos 2 (x) has an inversion symmetry about x = 
and we can always choose the instantaneous eigenbasis 
\v k {t)) to have real coefficients when expanded over plane 
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waves. As a result the above term goes to zero and the 
time derivative terms can be excluded. From Eq. (36) we 



can see that the last term in Eq. ( A3 ) , which depends on 



the light field fluctuation, can couple the different modes 
Sbj(t). At this juncture the coherent state approximation 
is made which assumes the light field fluctuations are 
independent of the atomic fluctuations i.e. 



8a(t) « d(t) 



2/i / dre- jA(t)(t ~ r) |(T). (A4) 
o 



This ensures that the different modes are uncoupled and 
we can write down their independent dynamics as 



dSbj 
~df 



(A5) 



where 



Uj(t) = iVNU (vj(t)\P{t) cos 2 (x)\<p(t)) (A6) 
P{t) = (a*(t)d(t) + a(t)$(t)) ■ (A7) 



These equations describe the atomic fluctuation dynam- 
ics in terms of a collection of independent oscillator 
modes that are acted upon by the shot noise force P(t). 
As described in [63] . we can now use Fermi's golden rule 
to derive a rate equation for each of the oscillator occu- 
pation numbers SNj(t) = (Sbj6bj) 



d(SNj 
dt 



= (T uj -T dj )(SN j )+T v 



(A8) 



where the damping and diffusion rates are 



These depend on the spectral density (power spectrum) 
of the shot noise force 



2nn 



(A 



off 



(A9) 



In the above expressions the shot noise spectrum is eval- 
uated at the shifted oscillator frequencies defined by 
Uj = Ej(t) — /i(t) with the instantaneous chemical po- 
tential fi(t) = {(p{t)\H{t)\(p{t)) . This shifting helps in 
removing the slow time dependence of the couplings Uj 
(derived from the meanfield Bloch oscillations) when de- 
veloping the rate equations. Since the total number of 
atomic excitations 8N(t) = ^2^SNj(t), Eq. (A8| gives 
an approximate description of the atomic fluctuation oc- 
cupation number. Since the damping and diffusion rates 
for the different oscillators are not the same, it is in gen- 
eral not possible to write down an equation similar in 
form to Eq. dA8) for SN(t). 



As mentioned in the main part of the paper in Section 



VII the system can reache a dynamic steady state in the 



sense that SNj is slaved to the instantaneous values of the 
transition rates. In this case, we can put d{SNj)/dt = 
and thus 



Weq(t) 



T dj (t)-T uj (t) 



4Af(t)^(i) 



(A10) 



where the negative sign is fine because AJ; ff (t) is negative 
on the cooling side of the effective resonance. 



Appendix B: Two time correlation calculation 

When the intracavity light field is written as d(t) 
a(t) + 5a(t), the signal variance is given by: 



(AA> 2 (cj,T)) = n 



cos 



\ujt) (\a{t)\ 2 + (6a<6a(t))) dt 



2k 1 Re 



dt\dt<x cos(wti) cos(u}t 2 )a(ti)a(t 2 )(8d ji (ti)5cr (i 2 )) 



dt\dt2 cos(a;£i) cos(wi2)a;* (ti)a(t2){5a(ti)Sa' (ta)) + / dtidt2a(ti)a* (£2) cos(a>ti) cos(u)t2)(5cr (ti)5a(t2)) 



2( K f/ 2 Re 



dtidt 2 a* (*i)a(t 2 ) cos(wti) cos(uit 2 )(^ b {ti)Sd' ! ' (t 2 ))9(t 2 - h) 



(Bl) 



In this appendix we provide details of how we numeri- 
cally compute the signal variance (and hence the SNR). 
The important extra computational step compared to the 



covariance matrix calculation in Eq. (31) is the evalua- 



tion of the two time correlations such as {8w (ti)5<v fa)) 
In the vector notation for the fluctuations, the two 
time correlations arc elements of the correlation matrix 



A(ti,i 2 ) = (^(ti) J R T (i 2 ))- The time evolution for the 
correlation matrix is given by: 

i^A(i,io) = M{t)A(t 1 t ) + i(Z(t)R T {t )). (B2) 
dt 

Let us consider the case when t > t . Then the last 
term in above equation gives a correlation between the 
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Langevin operators at some future time t and the system 
fluctuation operators at to. Due to the delta-correlated 
nature of the Langevin noise this term will be zero. This 
means that Eq. (B2| becomes homogeneous and we can 



solve it with the initial condition at t — to, A (to, to) — 
C(to). Also note that the time evolution operator for the 
numerical evolution in Eq. ( B2 ) is same as the one for the 



covariance matrix [denoted by G(t) in Eq. (31)], which 



is an expression of the quantum regression theorem [68] . 
A separate computation for t < to is not needed since 
they are related to the elements of A(i,t ) with t > t 
by complex conjugation. For example: 



(So! (t )6a(t)) = (<5a(i)5a T (t )>*. 

We can evaluate the correlation (£b(ii)<5a l '(i2)) using a 
similar approach as above for the time evolution of the 
vector (£&(£i)i?(t2)). In this case the initial condition for 
the evolution is (tb(to)5<$ (to)) — Since the evolu- 

tion operators for the correlation matrix and covariance 
matrix evolution are the same the calculation can be per- 
formed without additional computational cost. The main 



difficulty in computing the signal variance arises from the 
fact that the two time correlation functions are not sta- 
tionary. As a result, in order to evaluate the integrals in 
Eq. (Bl ) the correlation matrix needs to be computed for 



all values of < t\,t% < T. This is the memory intensive 
step in the computation and we simplify the situation 
by performing the correlation matrix computation over a 
coarser grid than the one used in the numerical solution 
of Eq. (31). This is justified since we find typically the 



correlation matrix elements do not change significantly 
over the very short time steps chosen in the solution of 
Eq. (31). Moreover, for the results presented in Section 
VIII[ we have taken care to check that the numerical so- 



lutions converge to a value independent of the size of the 
coarse grid. The necessity of evaluating two-time correla- 
tors over a two dimensional time grid is the main limiting 
factor to the maximum integration time for the SNR cal- 
culations. Another difficulty that we faced was that for (3 
values larger than the ones that we have presented here 
we found that the size of the coarse grid needs to be 
essentially matched with the size of the finer computa- 
tional grid over which Eq. (31 ) is solved. As a result the 



calculation becomes very memory intensive indeed. 
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